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Abstract 

Systems of interest in physics are usually composed by a very large number of interacting par- 
ticles. At equilibrium, these systems are described by stationary states of the many-body Hamil- 
tonian (at zero temperature, by the ground state). The reaction to perturbations, for example to a 
change of the external fields, is governed by the time-dependent many-body Schrodinger equation. 
Since it is typically very difficult to extract useful information from the Schrodinger equation, one 
of the main goals of non-equilibrium statistical mechanics is the derivation of effective evolution 
equations which can be used to predict the macroscopic behavior of the system. In these notes, 
we are going to consider systems of interacting bosons in the so called Gross-Pitaevskii regime, 
and we are going to show how coherent states and Bogoliubov transformations can be used to 
approximate the many body dynamics. 

1 Introduction 

In these notes, we review some recent results concerning an important question in non-equilibrium 
statistical mechanics which is the derivation of effective evolution equations approximating the dy- 
namics of many body quantum systems. 

When studying systems of interest in physics one typically has to choose between two different 
approaches. On the one hand, one can opt for a description based on the fundamental laws of 
physics. In this case, one would describe the system in terms of its elementary components and 
their interactions. So a chemical sample would be described as a system of many atoms, interacting 
through the laws of quantum mechanics. A galaxy could instead be described as a system of many 
stars interacting through Newtonian or Einstein gravitation. 

One can opt, alternatively, for an effective description of the system, which does not resolve the 
single components, but instead focus on their collective behavior. For example, the motion of a fluid 
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would then be described by the Euler or the Navier-Stokes equations, a superconductor through the 
Ginzburg-Landau theory. 

Of course the first approach is more precise, and typically relies on very few assumptions. The 
second approach, on the other hand, is less precise, it only has a small range of validity, but it is 
much more accessible to computations and focus exactly on those quantities which are measurable 
and of interest for the observers. 

The importance of developing simple effective theories, approximating the behavior of systems 
of interest has always been clear to physicists. In 1929, just after the development of quantum 
mechanics and of the Dirac equation, Dirac wrote: "The underlying physical laws necessary for 
the mathematical theory of a large part of physics and the whole of chemistry are thus completely 
known and the difficulty is only that the exact application of these laws leads to equations much 
too complicated to be soluble. It therefore becomes desirable that approximate practical methods of 
applying quantum mechanics should be developed, which can lead to an explanation of the main 
features of complex atomic systems without too much computation". 

Given their practical importance, a major challenge of mathematical physics is the rigorous 
derivation of effective equations starting from fundamental theories in appropriate limiting regimes. 
This should first of all justify the use of the effective equations, which are often introduced on the 
basis of heuristics or phenomenological arguments. Most importantly, a rigorous derivation should 
clarify the limits of applicability of the effective theories and should give bounds on the error of the 
approximation. 

We start our investigation from many body quantum mechanics. We consider systems of A'^ 
identical particles which can be described by a wave function Vat G L^iM.^^). We are going to 
assume the particles to obey bosonic statistics. We will assume, in other words, the wave function 
ipN to be symmetric with respect to any permutation of the particles, i.e. 

^Ar(x^l, . . . ,X^Ar) = iPn{xi, ■ ■ ■ ,Xn) 

for all permutation vr S S^. The wave function, normalized so that HV^A^Ib = 1, determines the 
probability density \ip]\fixi, . . . , xn)\'^ for finding particles close to xi, . . . ,xn- More generally, an 
arbitrary observable, given as a self-adjoint operator on the Hilbert space , can be interpreted 

as a random variable with probability law determined hy ipN through the spectral theorem. 
The dynamics of many body quantum systems is governed by the Schrodinger equation 

idttpN,t = HNtpN,t (1-1) 

for the evolution of the wave function ipN,t- On the r.h.s. of (II. ip . is a self-adjoint operator 
on L^(]R3^), known as the Hamilton operator (or simply the Hamiltonian) . We will restrict our 
attention to Hamilton operators with two-body interaction, having the form 

N N 
^N = Y1 (-^-. + ^ext(Xj)) + A ^ V{xi - Xj) . (1.2) 
j=l i<j 

The sum of the Laplace operators is the kinetic part of the Hamiltonian, generating the evolution of 
free particles. V^xt is an external potential and V describes the (two-body) interaction among the 
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particles (the parameter A G M is a coupling constant, introduced here for later convenience). The 
Schrodinger equation is a linear equation, and it always has a unique global solution in L'^{R^^), 
obtained by applying the one-parameter group of unitary transformations generated by to the 
initial wave function ipNfli i-6. ipN,t = e~*^'^*V'Af,o- What makes the Schrodinger equation interesting 
and challenging, is the fact that the number of particles N involved in the dynamics is typically very 
large; ranges from values of the order 10^ in extremely dilute samples of Bose-Einstein condensates, 
up to value of the order 10^^ in chemistry , up to even much larger values in astrophysics. For such 
values of A^, it is usually impossible to extract information from the many-body Schrodinger equation 
which go beyond its well-posedness. For this reason it is important to find effective equations which, 
on the one hand, are simple to solve, either analytically or numerically, and which, on the other 
hand, approximate the solution of the Schrodinger equation in the interesting regimes. 

A simple but non trivial regime, where it is possible to obtain an effective description of the dy- 
namics, is the so called mean field regime, which is characterized by the fact that particles experience 
a large number of very weak collisions. The mean field regime is realized by the Hamiltonian ()1.2p . 
when N ^ 1 and |A| <C 1, so that A^A remains fixed, of order one. This last condition makes sure 
that the total force acting on each particle, resulting from the many weak collisions, is of order one 
and therefore comparable with the inertia of the particle. Hence, to study the mean field regime, we 
will consider the dynamics generated by the Hamiltonian 

N N 

Hf = ^ (-A,.^, + VUxj)) + ^ 5] V{x, - X,) (1.3) 

in the limit of large A^. 

To understand the emergence of a simple effective equation approximating the many body dynam- 
ics in the mean field regime, consider the evolution of a factorized (or an approximately factorized) 
wave function 

N 

^pN,oixl,...,XN) ^ JJv'(xj) 
i=i 

for a if £ L^(M^). Since particles interact, factorization cannot be preserved by the time evolution. 
Nevertheless, because of the mean-field nature of the interaction, one can still expect factorization to 
be approximately restored in the limit of large A^. We can expect, in other words, that, in the limit 
of large A^, the solution of the Schrodinger equation ipN^t = e~^^N *^7v,o can still be approximated, 
in an appropriate sense, by 

N 

i=i 

Under this assumption, it is then simple to obtain a self-consistent equation for the evolution of (ft- 
In fact, (II. 4p implies that particles are approximately independent throughout the evolution. The 
law of large numbers then implies that the total potential experienced by the j-th particles can be 
approximated by the average of the potential w.r.t. the distribution \(pt\'^, i.e. 

1 ^ r 

- ^ Vixi - xj) c^Jdy V{xj - y)\My)\' = (V * \ipt\'){x,) . 
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Hence, if factorization is preserved in the limit of large N, the one-particle orbital (ft has to evolve 
according to the self-consistent nonlinear Hartree equation 



idm = (-A + Fcxt) ^t + {V* \iptf) ift. (1.5) 

In order to obtain precise mathematical statements about the relation between the many body 
Schrodinger evolution and the nonlinear Hartree dynamics, we need to clarify in which sense factor- 
ization is approximately preserved in the limit of large N. To this end, let me introduce the notion of 

(k) 

reduced densities. For k = 1, . . . , N , we define the /c-particle reduced density 7]^^ associated with the 



solution V'AT,* = e~*^™ *^7v,o of tlis Schrodinger equation by taking the partial trace of the orthogonal 
projection onto ipN,t over the last {N — k) particles. In other words, is defined as a non-negative 
trace class operator on Lp'iM?^) with integral kernel given by 



7^^{(xi, . . . , Xfc; , . . . , x'fc) = j dxk+i ■ ■ ■ cIxn iPnA^i^ • • • ' Xk,Xk+i, ■■■ , xn) 

^ V'Af,i(^'l) • • • ) ^fc! ■ ■ ■ ■> Xn) ■ 

(k) 

We use here the convention Tr 7]^ ^ = 1, for all = 1, . . . , A'^. For A; < A^, it is clear that the fc-particle 

(k) (k) 

reduced density 7)^ ^ does not contain the full information about the system. Still, knowledge of 7)^ ^ 
is sufficient to compute the expectation of any /c-particle observable, i.e. of any observable acting 
non-trivially on at most k particles. In fact, if O^^^^ is a self-adjoint operator on L^(R^'^) and if 
Qik) (g, -^(N-k) (|g]2otes the operator acting as O^'^^ on the first k particles and as the identity on the 
last (A^ — k) particles, we find 

^N,t, (0('=) ^ V^,) = Tr (0('=) = Tr O(^) 7i?,l. 

It turns out that the language of the reduced densities is the correct language to describe the 
convergence of the many body evolution towards the Hartree dynamics. 

Theorem 1.1. Let V be regular enough (see discussion below). Let tp^ = , for a £ H^(M.'^). 

Let '4)N,t = e~ be the evolution of ip^, as generated by the mean field Haniiltonian ^1.3\) . 

Then, for every A; G N, t S M, 

as N — 00, where ipt is the solution of the nonlinear Hartree equation il.5]) . with initial data 



The first proof of this theorem was obtained by Spohn [18] , for bounded interactions ||^||oo < 00. 
Spohn's approach was extended to potentials with a Coulomb singularity by Erdos and Yau [7j . 
These proofs were based on the direct analysis of the time-evolution of the reduced densities. From 
the Schrodinger equation for ^Nt, it is simple to derive a hierarchy of equations, commonly known 
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as the BBGKY hierarchy, describing the evolution of the reduced densities: 



+ 



N 



Eh 

Kj 



H 77— 2^ Trfe+1 V{xj - Xfc+i), 7]y_ 



(k+i) 
t 



:i.6) 



As — )• oo, the BBGKY hierarchy formally converges to the infinite hierarchy of equations 



K 



'A + ycxt(a;] 



(fc) 



+ ^Trfe+i 



(^•+1) 

' /oo,t 



:i.7) 



for all A; G N. It is then simple to check that this infinite hierarchy of equations has a factorized solu- 
tion = {kpi\®^ ^ given by products of solutions of the Hartree equation (jl.Sp . This observation 
suggests a strategy to obtain a proof of Theorem ll.il First, one has to show the compactness of the 
sequence of families {7^^}^^ with respect to an appropriate topology. Secondly, one characterizes 

the limit points {7^\}fe>i of the sequence of reduced densities as solutions of the infinite hierar- 
chy (ll.7p . Finally, one shows the uniqueness of the solution of the infinite hierarchy. This technique 
is very powerful, but it has few downsides. In particular, it does not give much information about the 
rate of the convergence towards the Hartree dynamics; this is an important point, if one is interested 
in applications to real physical systems, where the number of particles N is very large, but of course 
finite. Without a bound on the rate of the convergence, it is impossible to know whether N is large 
enough for the Hartree equation to be a good approximation. 

Motivated by these considerations, in a joint work with I. Rodnianski [T7] , we gave a new proof 
of the convergence towards the Hartree dynamics for potentials with Coulomb singularities. Our 
proof was based on ideas introduced in a work of Hepp [12] and then extended by Ginibre-Velo [9] 
. It had the advantage, compared with previous proofs, of giving precise bounds on the rate of the 
convergence (the technique was generalized in a joint work with L. Chen and J.-O. Lee [2] , to obtain 
optimal bounds). 

More recently, Knowles and Pickl [15] developed a different approach which also covered more 
singular potentials. Frohlich, Knowles and Schwarz |8j proposed an alternative point of view, express- 
ing the convergence towards the Hartree dynamics as a Egorov-type theorem. Grillakis, Machedon 
and Margetis [101 ITT] obtained a stronger convergence by considering a more general approximation 
of the full evolution in the mean field limit. 



2 The coherent states approach to mean field evolution 

In this section, we are going to briefly review the approach developed in collaboration with I. Rodni- 
anski [T7] . The starting point of the analysis are ideas introduced by Hepp [I2] to study the classical 
limit of quantum mechanics, which are based on the analysis of the evolution of coherent states in 
the Fock space. 
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The bosonic Fock space T over L^(R^) is defined as the direct sum 



where L^(M^") denotes the subspace of L^(]R'^") consisting of ah functions symmetric with respect to 
permutations of the n particles. Vectors in T arc hence sequences {V'^"^''}n>i, with V^''^ G L2(m3"). 
On T we can describe states where the number of particles is not fixed. We introduce the number 
of particles operator AA, which is defined by 



n>0 ■ 



Eigenvectors of H have the form {0, . . . , 0, ^^^^ , 0, . . . } with only one non- vanishing entry; these 
vectors describe states with fixed number of particles n. An important example of such a state is 
the vacuum vector $7 = {1, 0, . . . }, describing a state with no particles at all. 

It is very useful to introduce creation and annihilation operators. For / G L^(R3), we define the 
creation operator a*(/) and the annihilation operator a(/) by 

(a* (xi, . . . = ^ . . . . . . , 

(a(/)*)(")(xi,...,ar„) = V^y" V("+^)(x,ari, . . . 

for * = {V'*^"^}n>o £ ^- It is easy to check that a(/) is the adjoint of a*{f). Creation and annihilation 
operators satisfy canonical commutation relations 

[a{f),a*{g)] = {f^g)^^ [a{f),a{g)] = [a* {f) , a* {g)] = 0. 

We also introduce operator valued distributions a*,ax, defined so that 

«*(/) = j dxf{x)al and a{f) = J dxf{x)ax 



for all / G L^(R3). Expressed in terms of these operator valued distributions, the number of particles 
operator is given by 

J\f = J dxa^Qx ■ 

It is important to observe that, although creation and annihilation operators are unbounded, 
they can be bounded with respect to the square root of M. More precisely, we have 

||a(/)*|| < j dx \f{x)\ \\a,n < II/II2 (^j dx ||a,*f ^ ' = Wfh ||A^^/'*|| (2.8) 

and similarly 

||a*(/)*||<||/||2||(Ar+l)V2^r||. (2.9) 
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Next, we introduce a Hamilton operator on the Fock space J-. To simplify a bit the notation 
we will neglect the external potential and let 



Un = J Vxa%Vxax + j ^^'^V ~ ala^ayttx ■ (2.10) 

Since every term in T-L]^ has the same number of annihilation and creation operators, T-L^ commutes 
with the number of particles operator M. This implies that T-Ln leaves each sector with fixed number 
of particles invariant. It is simple to compute the action of Hn on these sectors. If = {ip^"'^}n>o, 
we have 



with 



j=l i<j 

In particular, we observe that, when restricted on the A^-particle sector, the Hamiltonian T-L^ coincide 
exactly with the mean field Hamiltonian (jl.3p (with 14xt = 0). This implies that the evolution of 
a Fock space vector ^' = {0, . . . , 0, ■0Af , 0, . . . } is exactly the mean field evolution discussed in the 
previous section, i.e. 

e-^^-*{0, . . . , 0, Vtv, 0, . . . } = {0, . . . , e-^5^?'V7V, 0, . . . } . 

The advantage of a Fock space representation is that we can consider more general initial states. In 
particular, we want to study the evolution of so called coherent states. 
For a if in L^(M^), we define the Weyl operator 

The coherent state with wave function ip is then described by the Fock space vector W{(p)Q. A 
simple computation shows that 



Hence coherent states do not have a fixed number of particles. Instead they are given by linear 
combinations of states with all possible number of particles. Since W*{ip) = W{—ip) = W~^{ip), 
Weyl operators are unitary, and therefore coherent states are always normalized ||VF(93)r2|| = 1. An 
important observation is that coherent states are eigenvectors of all annihilation operators. In fact, 
it is simple to check that 

W*{ip)axW{ip) = ax + ip{x), and W*{ip)alW{ip) = al +7p{x) (2.11) 

and therefore that 

W*{v)a{f)W{^) = a{f) + (/, and W* {ip)a* {f)W {^) = a*{f) + (99, /) 
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for every / e L^(M^). Since a{f)Q = for all / G L^(M^), these equations imply that 



aif)w{^)n = w{ip)w*{ip)aif)wi^)n = w{ip){a{f) + {f,^))n = {f,^)w{^)n 

confirming that W{(p)i^ is an eigenvector of a(/), with eigenvalue (/, 97), for any / € L^(M^). As we 
will see shortly, the algebraic properties of coherent states and Weyl operators substantially simplify 
the study of their dynamics. Using ()2.1ip . we can also compute the expectation of the number of 
particles in the coherent state W{(p)Q. In fact 

{w{<f)n,j\fw{ip)n) = j dx{n, « + (^(x))(a^ + ip{x))n) = WcpWl ■ (2-i2) 

Similar computations show that the number of particles in the coherent state W{ip)0, is a Poisson 
random variable with mean and variance \\^\\'^- 

Next, we study the time-evolution of initial coherent states, as generated by the Hamiltonian 
()2.10p . In order to recover the mean field limit discussed in the previous section, the number of 
particles must be related with the parameter appearing in (I2.10p . Although we cannot ask the 
initial coherent state to have exactly A^ particles, we can at least require the expectation of the 
number of particles to be A^. We fix therefore a ip £ L^(M^), with \\ip\\2 = 1, and we consider the 
time evolution 

of the initial coherent state W{^/N(p)Q. By (j2.12p . the expected number of particles is given by 
-^llv'lli ~ desired. In particular, we are interested in the reduced densities associated with '^N,t- 

To make the presentation simpler, let us focus on the one-particle reduced matrix F^^^ associated 
with "^N^f It turns out, that the integral kernel of is given by 

Since M is preserved by the time evolution, we have 
for alH G M. Hence 



^ {aye-'''^^^W{VN^)n, a^e-''^^*W{VN^)n) . (2.13) 



Because of the mean field character of the interaction, we may expect the evolution of the initial 
coherent state to be again approximately coherent, with an evolved wave function ipt, obtained by 
the solution of the Hartree equation (|1.5p . i.e. 
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If this is true, e~^'^'^^W{y/N(f)fl should be approximately an eigenstate of the annihilation oper- 
ators ax,ay, with eigenvalues ^/N ipt{x) , ^fN ^t{y) ■ Motivated by this observation, we expand the 
annihilation operators ax, ay on the r.h.s. of (I2.13p . around their mean field values. We find 

= '^t{x)Tpt{y) 

+ 4 Uay - VnMv)) e-'^'''W{^ip)n, {ax - VnMx)) e-'^^'WiVNip)n) 

N \ I (2.14) 



N 



N \ 

We observe now that the fluctuations {ax — ^/N(pt{x)) and {ay — ^/N(pt{y)) can be obtained by 
conjugating ax and ay with evolved Weyl operators, i.e. 



{ax - VN^t{x)) = W{VN^t)axW*{VN^t), 

and similarly for the fluctuations of ay. We define the fluctuation dynamics as the two-parameter 
group of unitary transformations 

UN{t;s) = W*{VN^t)e''^''^'-''^W{VN^,) . 

Then (|2.14p can be written in the compact form 

rSiJi(x;y) - ipt{x)lpt{y) 

= ^{n,U*^{t;0)alaxUN{t;0)n) 

+ ^^{n,u*^{t; 0) a*UN{t; o)n) + ^^{n,u*^{t; 0) axUN{t; o)n) . 

\/N ^ \/N 



The term ipt{x)(pt{y) on the l.h.s. is just the integral kernel of the orthogonal projection \ipt){ipt\- 
Hence, to show the convergence T^''^ — )• \ipt){ipt\, it is enough to bound the error terms on the r.h.s., 
and prove that they vanish in the limit of large A^. To this end, we recall from (j2.8p and ()2.9p that 
creation and annihilation operators can be bounded with respect to the square root of the number 
of particles operator M. The problem of proving the convergence — )■ \ipt){ipt\, in the trace norm 
topology and with an explicit bound on the rate of the convergence, reduces therefore to establishing 
bounds on the growth of the expectation 

{n,i(^{t-o)Mi(N{t;o)n) 

holding uniformly in N. 

To this end, we observe that the fluctuation dynamics U]\f{t; s) satisfies the Schrodinger equation 

idtUN{t]s) = CN{t)UN{t]s) 
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with the initial condition hlN{s; s) = 1 and with the time-dependent generator 
^N{t) = j dx V^^a* V^a^ + j dx {V * \'ft?'){x)al.ax 
+ / dxdyV{x -y)ipt{x)■ip^{y)alay 



+ ^ j dxdy V{x - y) {(pt{x)ft{y) ala*y + Lpt{x)Lpt{y)a^ay) (2.15) 
+ j dxdyV{x - y) a* {Tpt{y)ay + ipt{y)a*y) a^- 



+ ^y" dxdyV{x-y)alalaya^. 



In contrast with the original Hamiltonian (j2.10p . the generator Ci\i{t) contains terms where the 
number of creation operators does not match the number of annihilation operators. These terms 
do not commute with the number of particles operator M. As a consequence, the expectation of M 
is not preserved along the evolution Un- This is hardly surprising, since Un defines a fiuctuation 
dynamics, and fluctuations are expected to grow during the evolution. Although the expectation of 
the number of particles M is not constant, under appropriate assumptions on the regularity of the 
potential, which are for example satisfied for the Coulomb case V{x) = ±l/|x|, it was shown in the 
paper with I. Rodnianski [17j that, for every A: G N, there exists C,K > such that 

^,U*N{t;0){Ar + l)''U{t;0)^) <Ce^^'^{^,{Af + lf''+^^) (2.16) 



for all t S M, ^ G J-. The convergence towards the Hartree dynamics for the evolution of coherent 
states is then a simple corollary. 

Theorem 2.1 (Mean field evolution of coherent states). Suppose the potential V satisfies the operator 
inequality 

V^{x) < D{l-A) 

for a constant D > 0. Let ip e H'^{R^) and 'i!N,t = e"*^^*VF(\/iV(/?)i7. Let T^j^\ denote the reduced 
density matrix associated with \I' N,t ■ Then there exist C, K > such that 



Tr 



for all t G R. 
Remarks: 



Since (j2.16p bounds the growth of every power of the number of particles operator, we also 
find, similarly to (|2.17p . that for every A; G N there exist C, > such that 



Tr 

for all t G M. 



^f,-wt){^tr <^ (2.18) 
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From (|2.16p . we see that the convergence ()2.18p can be extended to the evolution of initial data 
of the form W{\/Nip)^, for arbitrary ^ G J" with {M + l)2'=+2\j/) < qo. 

Writing 



\\PNW{VNip)n\\ 

where P/v is the orthogonal projection onto the A^-particle sector, the bound (j2.16p can also be 
used to show the convergence towards the Hartree dynamics for factorized initial data. In this 
case, the analysis is a bit more complicated; optimal bounds on the rate of the convergence 
have been obtained in collaboration with L. Chen and J.-O. Lee [2j . The same techniques 
can be applied to initial A^-particle states of the form P]\[W{VN(p)^ /\\PnW{^/Nip)^\\, for 
arbitrary ^ £ with a{ip)^ = 0. 

It is worth observing that the coherent state approach discussed above not only implies the con- 
vergence towards the limiting Hartree evolution; instead it also give information about the behavior 
of the fluctuation dynamics in the limit of large A^. In fact, it was already proven by Hepp [T2] 
and Ginibre-Velo [9j that, as N ^ oo, Uiy{t;s) — )• Uoo{t;s) strongly, where the limiting fluctuation 
dynamics ZYoo(i; s) is defined by the Schrodinger equation 

idtUoo{t; s) = Coo{t)Uoo{t] s) 

with Uoo{s; s) = 1 and with the time dependent generator 

^oo(i) = j dx V^iC* V^a^ + j dx (V * \ipt\'^){x)alax 

+ / dxdy V{x - y)(pt{x)Tpt{y) a^ay 



+ dxdyV{x-y) [Lpt{x)>ft{y) ala*y + ipt{x)Lpt{y)a^ay) . 

The fact that £00 (i) is quadratic in creation and annihilation operators implies that the limit- 
ing fluctuation dynamics Uoo{t;s) acts as a Bogoliubov transformation. For f,g £ L^(M^), we let 
A{f,g) = a{f) + a*{g). Then we have 

A*{f,g)=A{J{f,g)), with J = ^ 

where J : L^(M^) — )■ L^(M'^) is the antilinear map defined by Jf = f. The canonical commutation 
relations take the form 

[^(/i,gi),yl*(/2,ff2)] = ((/i,5i),5(/2,52)W2 withS= r J 



A Bogoliubov transformation is a linear map 9 : ^^(M^) © L'^{M.^) -)> ^^(M^) © L'^{R^) with the 
properties 9J' = J9 and S = 9*S9. These conditions imply that, with B{f,g) := A{9{f,g)), we 
have 

B*{f,g)=B{J{f,g)) and 51), i?*(/2, 52)] = ((/i, ffi), 5(/2, 52)). 
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One can show the existence of a two-parameter group of Bogoliubov transformations 6{t; s) 
L2(m3) l2(m3) ^ l2(i^3) L-2(M.^)^ guch that 



W^it; s)A{f, g)Uoo{t; s) = Ai9{t; s){f, g)) . 
These maps satisfy the evolution equation 

idt9{t;s) = e{t;s)A{t) 
with 9{s; s) = 1 and with the time-dependent generator 

A{t) = 



Dt 
Bt 



where the hnear operators Dt, Bt : L' 



JBtJ 
JDtJ 



are defined by 



Dtf = -A/ + (V * \^t\^)f + (V * Tptf)^t 
Btf = iV *Tptf)-iPf 

Using this information about the hmiting dynamics Uoo{t]s), it is possible to prove a central 
limit theorem for the quantum fluctuations around the mean field evolution. Consider a factorized 
A^-particle initial data V'Af = 'f'^^ ^ for some (p G L'^iM?), and let il^N,t = eT^^N^Tp^ be its time 
evolution, as generated by the mean field Hamiltonian (|1.3p . with 14xt = for simplicity. Let O be 
a self-adjoint operator on L^(M'^), and denote by 0'^^\ for j = 1,. . . ,N, the operator on L'^{M.^^) 
acting as O on the j-th particle and as the identity on the other (N — 1) particles. At time t = 0, 
the observables 0^^\ . . . , O^^-* define a family of independent and identically distributed random 
variables. For t 7^ 0, on the other hand, 0^^\ . . . , O^^^ are no longer independent. The convergence 
(j2.17p easily implies that they still satisfy a law of large numbers. In a joint work with G. Ben Arous 
and K. Kirchpatrick [1] , we proved that, for every t E M, the random variables O'-^-', k = 1, . . . , N , 
also satisfy a central limit theorem: 



1 ^ 

(^OO) _ (^^, o^t)) ^ Gauss (0, a. 



in distribution, as A — )■ 00. So, the fluctuations around the mean field evolution are Gaussian, with 
a variance 



[Qt {J(pt,J(pt) ,@t {J(pt,J^t)) 



Qt {jLpt,Jipt) , -j= (if,^) 



which can be expressed in terms of the Bogoliubov transfomrmation 9{t; s) discussed above. 
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3 The Gross-Pitaevskii regime 



In this section, we are going to consider a different regime in which effective equations can be derived 
from many body quantum dynamics. The motivation comes here from the study of the time-evolution 
of initially trapped Bose Einstein condensates. In typical experimental settings, the condensate is 
initially trapped by strong magnetic field. After cooling the gas at very low temperatures, the 
traps are switched off and one observes the resulting evolution. The goal is to provide an effective 
description of the dynamics of the initially trapped condensate. 

On the microscopic level, trapped Bose Einstein condensates can be described as systems of N 
bosons with the Hamilton operator 

N N 

hT" = (-^-. + ^ext(xj)) + N^V{N{xi - X,)) (3.19) 

i=l i<j 

where V^xt describe the trapping potential while the interaction potential V is assumed to be short 
range and repulsive (meaning that V > 0). The interaction in (j3.19p scales with the number of 
particles, so that its scattering length is of the order 1/A^. Recall that the scattering length of the 
potential V is defined through the solution of the zero-energy scattering equation 

-A + ^V^f = (3.20) 

with the boundary condition f{x) — )• 1 as |x| — )• oo. Since V has short range, it is clear that for large 

\x\, 

f{x) ~ 1 - — 

for a constant oq > which is defined to be the scattering length of V. Equivalently, 

Mao = j V{x)f{x)dx. (3.21) 

A simple scaling argument implies that f{Nx) is the solution of the zero-energy scattering equation 
for the rescaled potential N'^V{Nx); in particular, this shows that the scattering length of N'^V{Nx) 
is given by ao/N. 

Lieb, Yngvason and Seiringer [14] proved that the ground state energy Ejsi ol the Hamiltonian 
(j3.19p is such that 

lim — — = min £Qp((f) 

N-^oo N ipeL2{R3):||<^||2 = l 

where we introduced the Gross-Pitaevskii energy functional 

^^GP(^) = Jdx [\Vip{x)\^ + Kxt(x)|<p(x)|2 + 47rao \ip{x)\^] ■ (3.22) 



Hence, in the leading order, the ground state energy depends on the interaction potential only 

jtra 
'at 



through its scattering length. Lieb and Seiringer [13] also showed that the ground state of Hf^"^ 
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exhibits complete Bose-Einstein condensation in the minimizer (/)qp of the Gross-Pitaevskii energy 
functional (|3.22|) . In other words, the one particle reduced density associated with the ground 
state of (j3.19p is so that 

as — )■ 00 (in the trace norm topology). Hence, in the ground state of (|3.19p almost all particles, 
up to a fraction vanishing in the limit of large A^, are in the one-particle state described by the 
minimizer (f)Qp of ()3.22p . 

Suppose now that the boson gas is prepared in the ground state of H^j^^^ (by cooling the gas to 
almost zero temperature) and that afterwards the magnetic traps are switched off. It turns out that 
the Gross-Pitaevskii theory can also be used to describe the resulting evolution of the condensate, 
as generated by the translation invariant Hamiltonian 

N N 

Hn = Y, + Yl N^V{N{xi - xj)) . (3.23) 

j=l i<j 

In fact, it was proven in a series of works with L. Erdos and H.-T. Yau [3l[ll[5l[6] that, for any family 
of A^-particle wave functions i^N ^ L'^i^'^^), with finite energy per particle 



{i^N,HNM < CN 
and exhibiting complete Bose-Einstein condensation 

for a £ L^(M^), the evolved A''-particle wave function ipN,t = e~*^'^*i/'Ar still exhibits complete 
condensation, in the sense that 

7Nl^\'Pt){^t\ (3.24) 

as N ^ 00, where the dynamics of the condensate wave function ipt is determined by the solution of 
the time-dependent Gross-Pitaevskii equation 

idtipt = -A(pt + 87rao\ipt\'^ipt (3.25) 

with initial data (pt=o = ^- Another approach to show (j3.24p has been recently proposed by Pickl 

m ■ 

The Hamiltonian (j3.23p can be written as a mean field Hamiltonian 

N , N 



N 

j = l Kj 



with an interaction N'^V{Nx). As A^ — t- 00, at least formally, 

N^V{Nx) bo5{x) where 60 = / dxVix) . 
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Inserting the limiting potential bo6{x) into the nonlinear Hartree equation approximating the many 
body dynamics in the mean field limit, we find an equation similar to (j3.25p . but with a different 
coupling constant in front of the nonlinearity (69 instead of Svrao). The reason for the failure of the 
mean field analogy is that the two limits, the mean field regime on the one hand, and the Gross- 
Pitaevskii regime on the other hand, are very different from the physical point of view. In the mean 
field regime, there is a large number of very weak interactions among the particles. In the regime 
described by H]\[, on the other hand, collisions are rare (particles only interact when they are very 
close to each others, at distances of the order and, at the same time, very strong. As a result, 

the many body evolution generated by Hj^f produces a singular correlation structure varying on the 
same 1/A^ length scale characterizing the interaction potential. Correlations, which are negligible in 
the mean field limit, play here a crucial role; in particular, they are responsible for the emergence of 
the scattering length in the Gross-Pitaevskii equation ()3.25p . 

To better explain this point, let us consider the first equation in the BBGKY hierarchy for the 
evolution of the one-particle reduced density associated with the solution of the Schrodinger equation 
'ipN,t- Expressed in terms of the integral kernels of 7^ j and this equation takes the form 

= (-A^ + A^,,)7}JJ(x;a;') (3_26) 
+ (A^ - 1) j dx2 {N^V{N{x - X2)) - iVV(A^(x' - X2))) j'^^lix, X2; x', X2) . 

If we assume that condensation is preserved by the time evolution, the kernels of the reduced densities 
7^''^ and 7^ ^ should be approximately factorized. However, for large but finite A^, the kernel of 7^^^ 
should also contain short scale correlations between particles one and two. Assuming that these 
correlations can be described by the solution of the zero energy scattering equation, we obtain the 
ansatz 

lNlixi;x[) = ^ptixi)Tpt{x[) , 

7^^^j(xi,X2;Xi,X2) = /(A^(xi - X2))f{N{x[ - x'2))ipt{xi)ipt{x2)^ti.x'i)Tpt{x2) . 

Inserting this ansatz in (|3.26p . we obtain a self consistent equation for ipt, given by 

idm = -^Vt + {{N - l)N^V{N.)f{N.) * \^t?) . 

As N 00, ({3:21]) imphes that (A^ - l)N'^V{N.)f{N.) Siraod, and we indeed obtain the correct 
Gross-Pitaevskii equation (I3.25p . The presence of the correlation structure makes the proof of the 
convergence (j3.24p substantially more difficult, compared with the mean field regime. Although the 
original strategy proposed by Spohn [H] could still be applied, the arguments had to be modified in 
several points. 

The difficulty of the analysis developed in the papers with L. Erdos and H.-T. Yau [31 HI [5] , 
together with the fact that it did not provide any control on the rate of the convergence, motivated 
us to look for a different approach and, in particular, to understand whether the coherent states 
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approach discussed in Section [2] could also be applied in the Gross-Pitaevskii regime. Following 
naively the arguments presented in Section [21 we define the Fock space Hamiltonian 

Un = J dxVxa%Vxax + \ j ^^'^V ^'^Vi^ix - y)) a*a*ayax ■ 

We would like to compare the evolution generated by Hn with the Gross-Pitaevskii dynamics (j3.25p . 
For technical reasons, it is more convenient to compare it with the modified Gross-Pitaevskii dynamics 
governed by the equation 

idm = -A^t + {N^V{N.)f{N.) * \^t\^) ■ (3.27) 

As — )• oo, it is then simple to show the convergence of the solution of (j3.27p towards the solution 
of the original Gross-Pitaevskii equation (|3.25p . with an error of the order for every fixed time. 

Proceeding as in the mean field case, the comparison of the many body evolution with (j3.27p 
leads to the fluctuation dynamics 

l(Nit;s) = VF*(\/iV^t)e-^^^(*~"%(\/iV^,) . 

As explained in Section O the problem of showing the convergence towards the effective evolution 
reduces to the problem of controlling the growth of the number of particles with respect to the 
evolution Un- Computing the generator of Un, we find the two terms 



CN{t)= idtW*iVNipt) WiVN^) + W*{VN^t)nNW{VNipt) . 
In (j2.15p . the linear (in the creation and annihilation operators) terms arising from 



cancelled exactly with idtW*{\' N(ft) TV(v Ac/j). Because of the factor f{N.) in the nonlinearity of 
(j3.27p . the cancellation is now not complete. We find, with uj = 1 — f, 

CN{t) = ^ J dx(al [N'^Lo{N.)V{N.) * \ipt\^] {x)ipt{x) 

+ ax [N^co{N.)V{N.) * \ipt\^] (x^x)) ^^'^^^ 
+ higher order terms in a,a* . 

Hence, the generator CN{t) contains a large term, of order -v/iV, which does not commute with the 
number of particles operator M. Because of this contribution, it is impossible to prove bounds similar 
to (|2.16p on the growth of the expectation of the number of particles holding uniformly in N. 

The reason for the failure of this naive approach is that we are effectively trying to approximate the 
many body evolution through an evolved coherent state, where correlations are completely absent. In 
other words, we are trying to bound fluctuations around the wrong correlation-free state. To control 
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the fluctuations, we have to improve our ansatz for the many body evolution, taking in particular 
into account the correlation structure produced by the many body evolution. 
To this end, we introduce the integral kernel 

kt{x,y) = -Nuj{N{x - y)) ipt{x) ipt{y) 

where u = 1 — f, and / is the solution of the zero energy scattering equation ()3.20p . For |x— y| ^ ^/N, 
we have Nu}{N(x — y)) ~ ao/|a; — y| and therefore 

kt{x,y) ~ - "° ipt{x)(pt{y) ■ 
\x - y\ 

Note that kt is the kernel of a Hilbert-Schmidt operator. Using kt, we define 

The unitary time-dependent operator T{t) acts on creation and annihilation operators as a Bogoli- 
ubov transformation. In fact, we find 

r(t)a(/)T*(t) = a(coshfc,(/)) + a*(sinhfc,(7)) 



where cosh^j , sinh^t : L (M ) —5- L (R ) are bounded operators, defined by the absolute convergent 
series 



The singularity of kt{x,y) at x ~ y disappears when we consider the powers {ktktY , n > 1. For this 
reason, for most purposes, we can approximate coshfcj(/) ~ /, sinhfcj(/) ~ kt{f), and therefore 

T* {t)aif)Tit) c^aif) + a* (kj) , and T* (t)a* (/)r(t) ^a*if) + a{hj) . (3.29) 

We use now the unitary operator T{t) to implement the correct correlation structure, improving 
our approximation for the many body evolution. We introduce the new fluctuation dynamics 

UN{t; s) = T*{t)W*{VNipt)e-'^^^'~'^ W{VNips)T{s) . 

Then Z^Ar(t; s) satisfles the Schrodinger equation 

idtUN{t;s) = CN{t)UN{t;s) 

with a new time-dependent generator 

CNit) = [idtT*{t)] T{t) + T*{t) \idtW*{VN^t)\ W{VN^t)T{t) 
+ T*{t)W*{VN^t)nNW{VN^t)T{t) . 
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The derivative of the Bogoliubov transformation T{t) is harmless (because it only acts on the solution 
(ft of the modified Gross-Pitaevskii equation, and not on the singular correlation function uj{N{x — 
y))). Let us focus on the last two terms on the r.h.s. of (|3.30p . As in (j3.28p . we find a contribution 
linear in the creation and annihilation operators, which is now conjugated with the unitary operator 
T{t): 

N j dxT*{t)a^T{t) [N^uj{N.)V{N.) * I^Jtp] (x) (ftix) + hermitian conjugate. (3.31) 

From the last term on the r.h.s. of (|3.3U|) . we also find a contribution cubic in creation and annihi- 
lation operators. It is given by 

^ ^'dxdy N^V{N{x - y))My) T* {t)alala^T{t) + h.c. 
1 



N 
1 



dxdy N^V{N{x - y))^t{y) T*{t)alalT{t)T*{t)a^T{t) + h.c. 
dxdy N^V{N{x - y))^t{y) 

a* + j dzkt{x,z)a^ {^*y^ j dwkt{y,w)ay^ T*{t)axT{t) + h.c. 

where we used the approximation (j3.29p . Some of the terms arising from the conjugation with T(t) 
are not normally ordered (they have a creation operator on the right of an annihilation operator). 
To control these contributions, we have to bring them back to normal order. Consider for example 
the term 

^ ^' dxdydz N'^ViNix - y))iptiy)h{x, z)a,a;Titya,Tit) + h.c. 
1 



dxdydz N^ViN{x - y))^t{y) kt{x, z)ala,T{ty a^TH) + h.c. 
+ dxdy N^V{N{x - y))ipt{y)kt(x, y)T{tTa.T{t) + h.c. 



'N 

where we used the canonical commutation relation [a^jfl^] = 6{z — y). Hence, normal ordering the 
cubic terms generates the new linear term 

1 



dxdy N'^V{N{x - y))^t{y) h{x, y)T{tyaxT{t) + h.c. 

N I dxT{tyaxT{t) [N^V{N.)uj{N.)*\^t\'^] {x)Tp^{x) +h.c. 



which cancels exactly the large linear contribution (j3.3ip . More cancellations involve the quadratic 
and the non-normally ordered quartic terms. Using all these cancellations, and controlling the other 
terms appearing in Ci\i{t), we obtain uniform control on the growth of the expectation of the number 
of particles operator with respect to the new fluctuation dynamic U^. As a consequence, we obtain 
the following theorem, establishing the convergence of the many body evolution towards the Gross- 
Pitaevskii equation, with a bound on the rate of the convergence. 
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Theorem 3.1. Consider the time evolution 

^1;^^^ = e-''^^^W{^ip)T{{))^P 

for a tp ^ J- with 



Then there exist C,ci,C2 > with 



Tr 



c 

~ \^t){vt\ < ^=exp(ciexp(c2|t|)) 



(3.32) 



for every t G ^ 
Remarks. 



From the convergence of the one-particle reduced density T^^^ towards a rank-one projection, 
one also obtains convergence of the higher order reduced densities, with rate N~^/^. For every 
fc > 1, there exist C, ci,C2 > with 



Tr 



N,t 



\v^t){^t\'^^ < CiV-i/4exp(ciexp(c2|t|)). 



• Theorem 13.11 continues to hold if the Hamiltonian contains an arbitrary external potential 14xt 
(assuming —A -|- Vcxt to be self-adjoint). 

• The bound (j3.32p deteriorates fast in time because it depends on high Sobolev norms of the 
solution ipt of the Gross-Pitaevskii equation ()3.25p . Assuming these norms to remain bounded 
in time, the estimate (|3.32p would only deteriorate exponentially in time. 
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